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Abstract 



Recently, Dymarsky, Kuperstein and Sonnenschein constructed an embedding of flavor D7- 
and anti-D7-branes in the Klebanov-Strassler geometry that breaks the supersymmetry of 
the background, yet is stable. In this article, we study in detail the spectrum of vector 
mesons in this new model of flavor chiral symmetry breaking and commence an analytical 
study of the scalar mesons in this setup. 
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1. Introduction 



The AdS/CFT correspondence [1,2] and its many non-conformal generalizations (often re- 
ferred to as gauge/gravity or gauge/string dualities; for a recent review, see [3], and refer- 
ences therein) have become promising and powerful tools to study non-perturbative physics 
of strongly coupled gauge theories in general and QCD in particular. There have been ap- 
plications in the field of heavy ion physics, where experiments at RHIC have revealed that 
the QCD plasma is strongly coupled and close to a perfect quantum liquid (for a recent 
discussion see [4]). The gauge/gravity correspondence can be applied and help elucidate the 
strongly coupled particle and hadronic physics of the quark gluon plasma [5,6,7]. Another 
subject where gauge/gravity techniques have recently been sucessfully applied, are strongly 
coupled condensed matter systems (cf. [8,9] and references therein), such as quantum critical 
systems and high temperature superconductors [10,11,12]. 

One of the best studied examples in the context of top-down 1 gauge/string models is the 
confining, non-conformal background with constant dilaton introduced by Klebanov and 
Strassler [18], which is dual to a 3 + 1-dimensional M = 1 supersymmetric field theory with 
gauge group SU(N + M) x SU(N). This field theory arises as the near-brane limit of N D3- 
branes and M wrapped D5-branes at the tip of the conifold. It exhibits logarithmic running 
of the gauge couplings which, along the RG flow, undergo a cascade of Seiberg dualities, 
where N repeatedly decreases by M units. In the far IR, the duality cascade must stop 
and the IR dynamics removes the conifold singularity via a gluino condensate that breaks 
the exact anomaly-free R-symmetry /L 2 m down to Z 2 . We are left with a warped deformed 
conifold whose deformation parameter e 2 is related to the expectation value of the gluino 
condensate. The gauge theory dual, in the far IR, is in the same universality class as Af = 1 
SU(M) SYM theory, albeit with extra massive adjoint scalars and fermions (for a review, 
see e.g., [27]). 

The Klebanov-Strassler model and its generalizations are particularly interesting phenomeno- 
logically, since they can be utilized in many ways. An important role is played by the 
-D7-branes in this context. Besides providing a holographic realization of QCD-like flavor 
physics of mesons, which will be the main focus of this paper, and baryons, it is also a cru- 
cial ingredient in string model building, i.e., compactifications with warped throats [19,20]. 
Furthermore, it seems conceivable that cosmological inflation consistent with present obser- 
vational data can be realized in such models [21,22,23]. The existence of meta-stable vacua 
has been investigated in [24]. Spectra of glueballs in the deformed conifold theory have been 
obtained in [25] and in [26] a finite temperature version of the Klebanov-Strassler model was 
used to study the physics of the quark-gluon plasma. 

From a theoretical point of view, much progress has been made over the past years, e.g., in 
connecting the Klebanov-Strassler background to the Gopakumar-Vafa open/closed string 
conjecture (see [27] for a review, and references therein). Moreover, its baryonic branch has 
been related to SU(3) structure backgrounds ([28], also see [29,30,31] for some recent devel- 
opments) and the resolved warped deformed conifold. Its mesonic branch has been studied 
in [32]. 

1 There arc also many phenomcnological models of QCD, inspired by holographic ideas. These bottom-up 
models are gencrically known as AdS/QCD, see, e.g., [13,14,15,16,17]. 
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Recently, Dymarsky, Kuperstein and Sonnenschein [33] were able to construct a non- super- 
symmetric embedding of D7- and D7-branes that constitutes a geometric realization of fla- 
vor chiral symmetry breaking in the Klebanov-Strassler model. This is based on the original 
ideas by Sakai and Sugimoto [34,35] and on similar results by Kuperstein and Sonnenschein 
[36] obtained previously for the Klebanov-Witten model [37]. In an early attempt [39] to 
incorporate a configuration of flavor D7--D7-branes into the Klebanov-Strassler model, that 
in fact predates the Sakai-Sugimoto model, the full solution could not be obtained due to 
complications arising from the non-trivial NS-NS flux in the Klebanov-Strassler background. 
These problems have been overcome in [33] by generalizing the results in [36]. 
In a recent paper [40], the authors studied phenomenological aspects of pions and vec- 
tor mesons in the Kuperstein-Sonnenschein holographic model of chiral symmetry breaking 
based on the Klebanov-Witten background. The numerical results for the form factors, par- 
ticularly for the pions, indicated that the meson physics in that model very well reproduces 
experimental data and QCD results. It is therefore quite natural to study phenomenological 
aspects of the Dymarsky-Kuperstein-Sonnenschein model, which is dual to a (large N c ) QCD- 
like theory in the IR, rather than a quiver gauge theory with gauge group SU(N) x SU(N), 
since it is based on the Klebanov-Strassler background. 

The focus of the present article is the investigation of vector and axial-vector mesons (in- 
cluding pions), scalar mesons and their respective spectra within the framework of the DKS 
.D7/.D7-brane embedding in the Klebanov-Strassler background. The paper is organized as 
follows: In order to be self-contained, we review the construction of Dymarsky-Kuperstein- 
Sonnenschein in section 2. Vector and scalar meson spectra are investigated in detail in 
section 3. While we find numerically robust results for the vector mesons, we did not find a 
numerically stable eigensystem for the scalar mesons, thus we only present analytical results 
on the corresponding Sturm-Liouville problem. We conclude with some remarks and an 
outlook in section 4. 

2. Review of the Dymarsky-Kuperstein- Sonnenschein model 

The Klebanov-Strassler background is a special case of the so-called imaginary self-dual 
(ISD) backgrounds in type IIB superstring theory with constant dilaton e* = g s . We will 
first review the strategy for constructing non-supersymmetric _D7-.D7-brane profiles for the 
general case of ISD backgrounds before specializing to the Klebanov-Strassler case. 

2.1. Non-BPS D7-brane embeddings in ISD backgrounds 

ISD condition. Let us briefly recall the definition of ISD backgrounds of type IIB super- 
string theory with constant dilaton e* = g s [41,42]. Such a ten dimensional background is a 
warped product 

ds 2 [w) = h~ l/2 dx^ + h 1/2 ds 2 {6) , (2.1) 

of flat Minkowski space-time M 3 ' 1 and an (internal 2 ) Calabi-Yau manifold Y 6 with metric 
ds 2 6 -p where the warp factor h may only depend on the coordinates on Y 6 . The self-dual 

2 In mathematics, Calabi-Yau manifolds are understood to be compact Kahler manifolds with a Ricci-flat 
metric. In physics, one often allows for non-compact examples, such as the (singular) conifold, which is a 
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R-R 5-form is given by 

F {5) = F {5) + B {2) A F (3) = (1 + *io)dC (4 ), (2.2) 

where F^ = dC( 2 ), F^ = dC( 4 ) and C( 4 ) = g^h^dx A dx 1 A dx 2 A dx 3 . The complex 
3-form flux G( 3 ) = F^ + yH(3) is imaginary self-dual w.r.t. the metric on Y 6 , 

G(3) = i*6 C( 3 ), (2.3) 

and vanishes in the Minkowski directions. It was shown in [41] that if G( 3 ) is of type (2, 1), 
then the background preserves M = 1 supersymmetry with four supercharges. 
Let us consider the R-R form C( 6 ). Its field strength is defined as 

F(7) = *10-F(3) - C(4) A i?( 3 ), (2.4) 

with *io-F(3) = /i" 1 *6 -F(3) A dx° A dx 1 A dx 2 A dx 3 . Note that the equation of motion for C( 2 ) 
implies that dF^ = 0. Using the explicit form of Cu\ given above, and observing that the 
ISD condition on G( 3 ) (2.3) implies 

^(3) = -97 1 H H( 3) , and * 6 F {3) = gJ 1 H (3h (2.5) 

it follows that F^ = [43]. Therefore it is consistent to set Ciq\ = for the rest of the 
discussion. 

D7-brane embeddings in ISD backgrounds. In this setup, D7-branes extend along 
the Minkowski directions and wrap a 4-cycle £4 inside the internal space Y & . Finding 
a classical D7-brane configuration entails solving the equations of motion for the scalars 
parametrizing the embedding and the field equations for the gauge fields on the worldvolume 
of the brane. Finding a solution is difficult and untractable in general; the situation is made 
even worse by the existence of a non-trivial NS-NS field = dB^ 2 ) in the background. A 
substantial simplification occurs if one is after an embedding preserving some supersymmetry. 
Scrutinizing the /t-symmetry of the theory, the authors of [44,45] showed that the problem 
reduces to finding a Euclidean D3-brane wrapping the 4-cycle £4. More specifically, in order 
to preserve the supersymmetries of the background, the D7-brane embedding and the gauge 
field living on its worldvolume have to satisfy the following conditions: 

• The 4-cycle £4 is holomorphic, i.e., £4 is the zero locus of an equation that can be 
written purely in terms of holomorphic variables 3 . 

• The gauge-invariant field strength J 7 = P[B/ 2 d + 2na'F is of type (1, 1) and anti- 
selfdual (ASD) for D7-branes (selfdual (SD) for D7-branes), 

T = — * 4 J 7 , ( T = + * 4 J 7 ) . 

cone over the base T 11 . It is possible, however, to compactify the (warped deformed) conifold, i.e., embed 
it as a "warped throat" into a compact Calabi-Yau manifold. 

3 See appendix A for a discussion of different coordinate systems used to describe the warped deformed 
conifold. 



4 



Unfortunately, supersymmetric embeddings of D7- and D7-branes can be shown to be rather 
trivial in the UV, namely the solutions exhibit only a single branch. Therefore they are not 
useful for modelling flavor chiral symmetry breaking. 
One of the main achievements of [33] was the following observation: 

Given any induced worldvolume metric g^, the corresponding D7-brane (D7-brane) ac- 
tion is bounded from below and minimized for J 7 anti-selfdual (selfdual). The resulting 
configuration is stable. 

Stated differently, the equations of motion for the D7-brane are solved by any embedding 
that extremizes the volume together with any ASD gauge field. The two questions are in- 
dependent of each other and can be addressed separately. This leads to the realization that 
not all ASD solutions are supersymmetric: 

If the embedding has to satisfy boundary conditions (e.g., invariance under a certain symme- 
try group) which are incompatible with holomorphicity but minimizes the volume J" s \/\g(4)\, 
for a certain induced within this class of geometries, then the resulting ASD solution 
is non-supersymmetric (sometimes also called non-BPS) and stable, since it minimizes the 
action. 

With this recipe in mind, we will now review the construction of non-BPS D7-D7 con- 
figurations in the Klebanov-Strassler background. 



2.2. Non-BPS D7-brane embeddings in the Klebanov-Strassler background 

The Klebanov-Strassler model. Klebanov and collaborators successfully generalized 
the theory of N D3-branes sitting at the tip of a singular conifold (the Klebanov- Witten 
model [37]) to include M fractional D3-branes (D5 branes wrapped on the collapsed 2-cycle 
of the conifold). In a sequence of articles [46,47], culminating in [18], they developed an 
understanding of the supergravity solution and its dual gauge theory. The ten dimensional 
supergravity solution is a warped deformed conifold (2.1), where h solely depends on the 
dimensionless radial coordinate r. The six dimensional metric on the deformed conifold is 
given by [28] (see also [55]), 



e 



4/3 



ds( 6 ) = — —K(t) cosh r 
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4 + e\ + e\ + e\ + — r-(ei£i + e 2 e 2 ) + " — (dr 2 + e 3 

coshr SK (t) 6 coshr 



2 

— — (dr' + 
coshr 

(2.6) 

where the vielbeins e iy i = 1,2 on S 2 and Ej, j = 1,2,3 on S 3 are defined in appendix A. 
Note that £3 = £3 + cos#id0i. Moreover, 

, , (sinhr coshr — r) 1 / 3 

K{T) = ^Tr ■ < 27 > 

The warp factor can be obtained by solving a six- dimensional Laplace equation, yielding 

h(r) = (^Ma') 2 2 2/ V 8 / 3 [°° dx ^otho;-l ^ _ 1/3 

J T sinh x 
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The 775-branes wrapped over the S 2 of the T 1 ' 1 are sources for magnetic R-R 3-form flux 
through the S 3 of T 1,1 . Thus, in addition to N units of 5-form flux, the solution also has M 
units of 3-form flux, i.e., 



F (3) = M, / F {5) = N. (2.9) 

5 3 JT 1 ' 1 



In the following the NS-NS Kalb-Ramond field 

5(2) = 6(r 

with 



1 / * 
e 1 A e 2 + e 1 A e 2 -\ : — (ei A e 2 + £i A e 2 , 

cosnr 



(2.10) 



6(t) = —-g s Ma V VV,, a _ -, (2.11) 



1 , , , cosh r(r cosh r — sinhr) 

-o s Ma' i 5— 

2 y 2sinh 2 r 

will play a crucial role. 

The dual gauge theory is a M = 1 £[/ (iV + M) x SU(N) gauge theory with two chiral 
superfields A 12 in the bifundamental (N + M, N) color representation and two chiral su- 
perfields 77i >2 in the bifundamental (N + M, N) color representation whose superpotential 
reads W = \e ik e jl AiBjA k Bi. The theory is invariant under a 377(2) x SU{2) x Z7(l) global 
symmetry and a discrete Z 2 m R-symmetry. The relative gauge coupling 

9T 2 ~ 9 2 2 ~ ^J 1 5(2) - , J s B {2) ~ M^ s ln(r/r uv ), (2.12) 

runs logarithmically due to the dependence of 77( 2 ) on the radial coordinate. Here, ruv is 
an arbitrary UV scale, where it should be noted that r 2 ~ p 4 / 3 for large radial distances, 
with p 2 = e 2 coshr (see appendix A). Thus, the magnetic 3-form flux is responsible for the 
breaking of conformal symmetry in this model. The self-dual 5-form flux becomes 

F {5 ) = (N + ag s M 2 \n(r/r lJV ))volT 1 '\ (2.13) 

where volT 1 ' 1 is the volume form of the conifold base T 1 ' 1 and a is a constant of order one. 
A discussion of the Seiberg duality cascade, both from the points of view of the gauge theory 
and the supergravity can be found in the original article [18] (cf. also [29]). There is an RG 
cascade, in which the 5-form flux F( 5 ), present at some UV scale Tuv, repeatedly drops by 
M units as j s2 77( 2 ) goes through a period, ultimately vanishing in the IR. The important 
observation is that this cascade can be related on the dual gauge theory side to a cascade of 
Af = 1 transformations on the gauge group factors, so-called Seiberg dualities. Since g\ and 
g 2 flow in opposite directions, these duality transformations have to be performed at scales 
Aj where one of the coupling constants diverges. Obeying certain matching conditions, it 
turns out that the new theory has the same form as the old theory under N — > N — M, 
thus resulting in a self-similar RG flow. When the cascade stops in the far IR, there are 
no 773-branes left other than the M fractional 773-branes, and the gauge theory is M = 1 
SU(M) SYM. The conifold singularity is removed by blowing up the S 3 of T 1 ' 1 , as can be 
seen by studying the resulting geometry via introducing a probe 773-brane: The probe brane 
lives on a deformed conifold. 

Non-supersymmetric D7-brane embeddings. We can now apply the strategy outlined 
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for the general case above to the case at hand. Recall that the D7-branes extend in the flat 
Minkowski directions. The 4-cycle S 4 describing the embedding on the deformed conifold 
necessarily extends along the radial direction r. Moreover, in order to be SU(2)r invariant, 
£4 must cover the S 3 of the T ,l base completely. Therefore, the D7-branes are located 
at points in S 2 that depend on r, similar to the situation in [36]. Hence, the embedding 
describes a trajectory (<K(t), #i(t)) on S 2 . Without loss of generality, one can set 61 = tt/2. 
The resulting induced metric = P[gn \] is the pullback of the ten-dimensional metric to 
the world volume of the D7-brane and reads 



ds 2 {8) = h- l ' 2 {r)dx^ + /* 1/2 (r)d4 4 , (2.14) 



where 



e -4/3 ds | = fM^ coshr 



4 



2 



tanh 2 r ( ^fi + — dr 2 + e 2 + ^ + " — e\ 

4 \\ \dr J 3if(r) 3 coshr / 1 2 3K (r) 3 cosh r 1 

(2-15) 

with 82 := £2 + 9 coJh t The classical configuration that solves the equation of motion can 
be determined from the DBI action of the D7-branes without gauge fields, integrated over 
the Minkowski and S 3 directions, 

1 /2 

Mo) |27r 2 8/ . rl f°° , coshr / 3A'(r) 3 sinh 2 r /<90i\ 2 \ 

SZ = W jj^VM-*. I (l + ^ (■£ ) J (2.16) 

yielding 

dcpi \ 2 coshr / K(t) sinh 2 r coshr \ 

9r"J " 3A(r) 3 sinh 2 r V A(r ) sinh 2 r cosh r / ^ ' ' 

This leads to the U-shaped configuration that is needed for a geometric realization of flavor 
chiral symmetry breaking: The D7-brane profile stretches down to a minimal radius r = tq, 
where jumps from —00 to +00, i.e., the D7-brane turns around at r = r , becoming 
a D7-brane. The minimal radial distance from the tip of the deformed conifold, tq > 0, 
parametrizes a one-parameter family of solutions to the embedding equations. 
We can calculate the total angle A0i covered by the trajectory on S 2 by integrating 

1 _ /2\ 1/2 r°° ^ ( coshr \ 1/2 / K(t) sinh 2 r coshr 

2 W Xo T V^(r) 3 sinh 2 ry V^(r ) sinh 2 r cosh r 



There are two interesting limits: For small r, K{j) approaches (2/3) 1//3 , thus in the limit 
t — > 0, the total angle evaluates to 

1 f°° dr 7T 

-A0 1 (r o ^O) = y — r-r^=2- (2-19) 



r 



This corresponds to an antipodal configuration. 

On the other hand, if r is large, K{r) « 2 1 / 3 e _T / 3 , and thus in the large ro limit, one obtains 

1 /6 

-A0 1 (r o ^cx)) = ^-7r, (2.20) 
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which exactly reproduces the result found in [36] in the conformal Klebanov-Witten setup. 
ASD gauge field on the D7-branes. To complete the construction of a non-supersymmetric 
.D7/-D7-brane embedding, we need to find an (A)SD gauge field configuration J 7 on its world 
volume. We will briefly reproduce here the results found in [33]. One starts with writing 
down a gauge-invariant combination of the field strength on E 4 in the A T = gauge, subject 
to the Bianchi identity, 

AF=P[ff (3) ], (2.21) 

with gauge field components on the S 3 , A a=5t& j. The pullback of the B( 2 ) field to the 4-cycle 
£4 is given by P[B/ 2 \] = b(r)ei A e 2 . Then the ansatz for the field strength reads 



Mr) u, ,\ , . ( , Mr) 



J 7 = P[B(2)} 

+A' 6 dr Ae 2 - (A 7 (t) - 6(r))ei A e 2 - A 5 (r)e 2 Ae 3 + A 6 (r)ei A e 3 , (2.22) 

where we have absorbed a dimensionful factor of 2ira' into the gauge field components. 
Here, the prime ' denotes derivation w.r.t. r. One should keep in mind that r is only a good 
coordinate on each branch separately. Using the metric (2.6), we can write down the ASD 
condition for the D7-brane branch, 

Mr) ~ = K(r)A 5 (r), (2.23a) 

A' 6 (r) = L T0 (r)A e (r), (2.23b) 

with 

T M ft 1 3 ^( r ) 3sinh2r ^^^ 1 2 A A(ro)sinh 2 rocoshr o y 1/2 
L - (T)= l 1+ 2 cosh r 0l(T) J = (, A(r)sinh 2 rcoshr J " ^ 

On the other branch, for the SD gauge fields on the Z)7-branes, one has to replace L T0 (r) — > 
—L to {t). It is also possible to switch to the globally well-defined coordinate 0(r) instead 
and reexpress the (A)SD equations in terms of <fi. It was shown in [33] that Aq = to ensure 
regularity of the solution (otherwise A§ oc e r for large r). Unfortunately, there is no analytic 
solution for the remaining two gauge field components. 

For r = 0, one finds that the only non-divergent ASD solution is given by A 5 (t) = and 
where Co is an integration constant and 
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It is necessary to set Co = to avoid an exponentially growing contribution to A-j{t) for large 
r. The non-divergent SD solution A 7 (t) on the D7-brane is 



A 7 (t) 



-S(t) 



26(1 



T 



df 

3K(f) 3 coshf 



Cl 



(2.27) 



where c\ ~ —0.83 has to be adjusted such that the SD and ASD solutions are continuous at 
the origin r = 0. We show the numerical solution for the background gauge field in figure 
1. For t > tq, it can be shown that there is a unique solution which is regular both in the 



A 7 (T) 




-10 



Figure 1: The background gauge field for r = 0: The ASD solution A-j{t) on the Z)7-brane 
branch (blue), and the SD solution A 7 (t) on the D7-brane branch (red). The constant c\ is 
chosen such that the two solutions match smoothly at the origin. 



UV and at the tip of the profile r = r (see appendix B of [33] for details). Note that the 
solutions are uniquely fixed by the boundary conditions at r = r , namely to ensure that the 
ASD gauge fields on the Z)7-brane branch are smoothly connected to the SD gauge fields on 
the D7-brane branch. 

3. Mesons 



3.1. Vector mesons 

According to the gauge/gravity correpondence, vector and axial- vector mesons correspond 
to fluctuations of the U(Nf) gauge fields living on the D7/D7-branes. The D7-brane action 
consists of two parts, 



S 



D7 



Sdbi + Swz 
-jj.7 I d 4 xdrd 3 {le~ l 



W) 



A e 



T 



(3.1) 
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where the first part is the Dirac-Born-Infeld action of the D7/D7 branes and the second part 
is the Wess-Zumino action. The induced metric g^) is given by (2.14) and F was defined 
above. For the case at hand, only the C(4) contributes to the Wess-Zumino action, 

Swz = -/i 7 J ^ A J 7 A (7(4). (3.2) 
Under the weak gauge field assumption, we can expand the DBI action as 

S DB1 [g,F + 5F] = SSIA9} + Sym[(5T) 2 } + 0[(5F)% (3.3) 

where ~ J^, -\/l^(4)T was use d above to determine the classical _D7-brane embedding 
(2.17). We want to expand this action around the classical action (3.1) and study a per- 
turbation of the form 5F = 2ira'F, where the only non-vanishing components of the field 
strength F are along the Minkowski directions and the radial direction r. Following [33], we 
define 



E:=g (8) +F yW\ = V 1^(4)1 + /r 1 V m (3-4) 
where ga\ is the induced metric on £4. Its inverse is (schematically) given by 

^ = f ftl/v "< 7i=£#Uf (%; - "- i/2 %^%0 1 (3.5) 



y/\9w \ + h 1 VW\ 

The action 5Sym ~ 0(F 2 ) can then be integrated over the S 3 directions to give rise to the 
five- dimensional effective action for the gauge fields (for details, see appendix C of [33]). 



St; >eS = -kJ ^xdr^\E\Xr{{E^F) 

= -kJ d 4 xdr (^C T0 (r)r)^ p F Xp F^ + D T0 (r)r,^F^A 
where k = 2{-Ka') 2 iijg^Yol S 3 . Here we have introduced 



(3.6) 



CVo(r) 



^ Hj^i cosh r + Mr) ^ (t) _ b{T) m^ + {Mr) _ b{T)) 4(T)> 

n A [ rr e 4/3 A» cosh r 

r (r) = y 1^(4)1^(4) = -r 7 (r) • (3-8) 



Since C to (t) depends on F, it is different for D7- and .D7-branes. Again, the two branches 
have to be smoothly connected at r = t . For the study of meson spectra below, it turns 
out to be advantageous to use rescaled, dimensionless variants of the functions introduced 
above to explicitly introduce a mass scale 4 M*: 



(3.9) 



4 As explained in [18], the four-dimensional mass scale is set by M g b ~ M* ~ g e ^ Ia , ■ This is the mass 
scale associated with the glueballs of the gauge theory. There is another independent mass scale associated 
with the mesons M meson ~ A^ g b/( T "o)j where /(r ) will be determined below. For r = 0, the two mass scales 
are identical, M mcson = M g b- 
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where k = {9aM 2 *' ? M 2 = 6 (g ^, )2 and 
a (r) = Hr)^M cosh r 

AT 



+A(r) ^(r) - S(r)^:J + (A 7 (t) - S(r)) 4(r), (3.10) 
n ^ /<(r) coshr 

Here we have defined the dimensionless functions 5 



oo 



h(r) = (g s Ma')- 2 e s/3 h(r) = 2 2 / 3 / dx XC ° th ^ 1 (sinh 2x - 2x) 1/3 , (3.12a) 

7 T sinh x 

u \ 2v ^ u \ /tt cosh r(r coshr- sinh r) ,„ 1I)U 

6 r = U > h r = ^"Ta ' 3.12b 

g s Ma' 2smh r 

= -^Wr). (3.12c) 

From now on we will drop the hats, but it is to be understood that we will work solely 
with dimensionless quantities. In the following, we want to study fluctuations of the gauge 
fields around the classical solution that correspond to vector mesons, axial-vector mesons 
and pions. For convenience, we will work in the A T = gauge. Then, can be expanded 

as 

oo oo 

A,(x,r) = V,(x) + Mx)Mr) + J2 v ^( x )^-i(r) + ^ ffl W(x)fe(r), (3.13) 

n=l n=l 

where the zero-mode ipo{r) will be addressed later and the fields and cffi are associated 
with vector and axial- vector mesons, respectively. In the A T = gauge the pion field U(x) 
appears in the expansion of 



V M (x) 



: ^(4>) + 47>)) (3-14) 

= ^+ (a l ,(x) + eg + \i- (a^(x) + eg c\ 

Mx) = \ (A%{x) - A%(x)) (3.15) 
= (A Lfl (x) + d») e; 1 - it. (A^(s) + <9 M ) C\ 

where A L/J (x) and A fiAt (x) are external gauge fields and we have introduced £jl = £_ := 

. n(z) 

e* At . The normalization conditions and equations of motion for the wave functions i^ n (T) 



5 Due to the different normalization of b(r) compared to [33], we will have a slightly different asymptotics 
for A-j(t), namely A^ir — s- oo) = \\ §(— 2r — 1), which differs by a factor of 6 -1 / 2 from the DKS result. 
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are determined to be 



k J dr Cr^r^mirjipnir) = 5m,n, (3.16) 
-{C T0 {r))- l d T {D T0 {T)d T ^ n {r)) = XnMr), (3.17) 

with X n := M 2 n /Ml 

Four-dimensional effective action for the mesons. Substituting the gauge field (3.13) 
into the action S^es, eq. (3.9), we find the four- dimensional effective Lagrangian for the 
vector and axial-vector mesons (disregarding divergent terms from non-renormalizable con- 
tributions), C^cs = £>4d,div + Y1T=2 £-4defi- The kinetic part of the Lagrangian reads 



41s = ^tr(0,# - d u v^f + ^tr(^) - d v ~a^f + tr(z^n + /^,) 2 

+M> (vM - |gVf + M> (a™ - (3-18) 



where, in order to diagonalize the kinetic terms, we have rewritten the meson fields in the 
following way: 

f'i-W' + lg:^ aW=.W + |^4 11 . (3.19) 

Moreover, we define 

M 2 n = X 2 n-xMl M 2 an = \ 2n Ml (3.20) 
V M = \{A L , + A Rll \ A,= l -{A L ,-A R ,). (3.21) 

The coupling constants between a massive vector meson v}? (axial- vector meson al ) and 
an external U(l) field V M representing a photon (an external axial U(l) field A^) are 

g v n = nM 2 n / _dTC T0 {T)tlj 2n -i{r), 

JD7/D7 

= -2kM 2 (D T0 (r)d T i; 2n ^(T)) , (3.22) 

r— >oo 

g an = nM 2 an [ _drC T0 (r)ij 2n (T)Mr)- (3.23) 

JD7/D7 

In the next section, we will proceed to numerically calculate the wave functions ?/Vi(t) in 
order to be able to compute the mass spectrum M 2 n and M 2 „. 

Study of the numerical wavef unctions. For r = 0, the (extremal) -D7/.D7-brane em- 
bedding touches the tip of the deformed conifold, where the S 2 shrinks to zero size. On the 
two branches the value of 0i is constant, i.e., 0' x = 0. At r = 0, <p[ diverges and the value 
of 0i jumps from — ~ to +|. In this case it is possible to describe the A 7 component of 
the background gauge field on the L>7-brane branch by a simple continuation of (2.25) to 
negative values of r (see [33]). Following [33], it is possible to introduce new coordinates 
similar to the ones introduced for the singular conifold [36]: 

y = t cos 0i r °=° 0, z = t sin 0i r °=° ±r, (3.24) 
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where the + sign corresponds to the D7-brane branch and the — sign corresponds to the 
D7-branch. It is clear that in order to study realistic particle physics scenarios, one needs 
to have tq > (e.g., to separate the mass scale of the mesons from the mass scale of the 
glueballs). However, we have not been able to obtain a numerical solution of the coupled 
differential equations for the background gauge field for tq ^ 0. Therefore, we will commence 
our investigation of the mass spectra for (axial-) vector mesons with the antipodal embedding 
corresponding to r = and leave the more involved case of general To as a future problem 6 . 
The general strategy is to employ a shooting method as discussed previously in, e.g., [39,40]. 
Antipodal configuration. In this extremal case, the functions Cq{t) and Dq(t) read 

Cb(r) = h(T)^coshT + (A 7 (r)-b(r))A' 7 (r), (3.25) 

A)(t) = ^4 cosh r, (3.26) 

where A 7 (t) is the rescaled version of the function given in (2.25). Moreover, L (t) = 1 and 
Lq(t) = 0. We can rewrite the equations of motion (3.17) as 

€(r) + ^rk(r) + ^KMt) = 0, (3.27) 

where 

D' (t) 3(rcothr - 1) - sinh 2 r 



Dq(t) 3(sinh r cosh r — r) 

Co(r) h{r) 2 (A 7 (r) - b{r)f 

D (t) K 2 (t) 3 K{rY cosh 2 r 



tanhr, (3.28) 

(3.29) 



We are now interested in finding the asymptotic behavior of the wave functions ip n ( T )- To 
this end, we expand them into a Frobenius series 

oo 

Mr) = ±e" QT ®n, m e^ mT . (3.30) 

m=0 

Note that 

"TdM = I < 331 » 

T->00 L)q(t) 3 

lim -^44 = 0. (3.32) 



A)(r) 



6 As discussed in some detail in [38], the To ~ or r$ ~ r e regime is of great phcnomcnological interest, 
since all meson masses are of the same order m g b and thus an approximate cancelation of the attractive 
and repulsive nuclear forces (as observed in nature) can possibly occur in this regime. A related important 
issue is the question whether the pseudo-Goldstone boson ct, which is parametrically lighter than the lightest 
vector meson for To >> 0, remains lighter in the To ~ regime, resulting in a net attractive nuclear potential 
at large distances. 
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From this behavior and imposing the normalization condition (3.16), we find that asymp- 
totically the normalizable solution must have a = |, i.e., it goes like 

^ n (r) = ±e-^ n (r), (3.33) 

where ip n (r — > oo) = const, at the boundary. It is easy to check that ^n( r ) satisfies the 
differential equation 

g(r) + X(T)i/ n (T) + (Y(t) + X n Z{r))i, n {r) = 0, (3.34) 

with 

X(r) = (3.35) 
D (t) 3 

F(r) = 9-3^)' (3 ' 36) 

Z(r) = gg. (3.37) 

Large r regime. For large r, we can expand X(r),F(r) and Z(r) as 

00 OO CO 

X(r) = ^(r)^, Y(T) = Y,yk(T)e~ 2k \ Z(r) = e-^^^(r)e-^, (3.38) 

fc=0 fc=0 k=0 

yielding 

2 16 8 , , 4 16 32 2 . . 

Xo[T) = --, xi(r) = -— + -r, x 2 (r) = -- — T + —t , . . . (3.39) 

/ x n / x 32 16 8 32 64 2 , o ^. 

Ito(r) = 0, yi (r) = T - -r, y 2 (r) = -- + -r - -r\ . . . (3.40) 

The expansion of Z{r) is a little more involved: Note that the expansion of k\t) OIU y 

contributes to the coefficients in the expansion, while the expansion of f k^Ios^ t P r °duces 
contributions to all the other terms as well. Then, 

, . 9/ o /15 3 \ . . 9/ o / 300 24 

2/3 /_23661 _ 203 31 A 

v 7 V 4 ooo ioo 5 y v ; 

By expanding 

oo 

Mr) = ^a„, fc (r)e^^, (3.42) 

fc=0 
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and using (3.34), one finds the following recursion relation for the coefficients a n> k'- 



k , k 2 



n,k + 



3 36 

m=0 



_ x 

S'rn I Ot n k— 12m ~~ m^n,k—V2m 



6 



+ A n ^ z m a n , k -3m-4: = 0. (3.43) 



m=0 

Setting a„ 5 o = 1, we obtain the first few coefficients as 

27 

«n,l = «n,2 = «n,3 = 0, a„, 4 = - — 2 2/3 A n (7 + 2t), 



729 

7 = 0, a n , 8 = X n (769 + 196r + 16r 2 ), 

3 r nl ~ 2 / 37665 3159 243 ,\ 
ct n 9 = a n in = a n ii = 0, a n \o = 1 h 2 1 X„\ 1 r H r 

n,9 n,W n,ll , n,12 g T g T „ ^ T ^ 2 56 / 



We can now repeat the same calculation for the D7-brane branch by using the SD background 
gauge field Aj(t) instead of Aj{t) and replacing the expressions for Cq(t) and Z(t) with the 
appropriate D7 expressions, namely 

7U\ Co h(r) 2 (A 7 (r) - Kr)f 
The appropriate large r expansion for Z(r) turns out to be 

oo 

Z(r) = e -^^z fc (r)e-^. (3.45) 



fc=0 



Again, we can employ a similar ansatz as before, ip n = ±e 3T ^> n with ip n = a n,ke 
and a n0 = 1. The corresponding recursion relation is given by 



(i 



k 

_T2J 



ii k i k ^ ^ f i 1 \ 

a n,k - ^ a n,k + + Xm ( a n,fc-12m — g(« — 12m)a ni fc_i 2m I + V m Oln,k-V2m 

m=0 L ^ ' - 

+ A n ^ = 0. (3.46) 

m=0 

Strategy. All this can now be used as input data to solve the equations of motion (3.17) 
numerically via a shooting technique. Since the equations of motion for the (axial-) vector 
mesons depend on the background gauge field A 7 (r) ^ A 7 (t), the coefficient functions are 
different for the D7- and D7-brane branches (cf. eqs. (3. 28), (3. 44)). Therefore we do not 
expect the wave functions ip n do exhibit parity symmetry about the origin, as was the case, 
e.g., in the Kuperstein-Sonnenschein [40] and Sakai-Sugimoto [34] models. However, in the 
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large r limit, the functions and j°ay , involving the background gauge field A 7 (t) / A 7 (r), 
tend to zero. It is therefore to be expected that the parity symmetry will be restored in this 
limit. As in e.g. [40], we are dealing with differential equations of the Sturm-Liouville type. 
We find it convenient to work with the z = ±r coordinate introduced above and denote the 
wave functions defined on the full domain of z by ty n (z). Then we will treat the problem 
of finding a discrete spectrum of (axial-) vector mesons living on the full D7/D7-brane 
embedding as a single differential equation on the complete domain of z, z e (—00, +00), as 
the natural extension of eq. (3.17), 

KM + wf\<(z) + An§^ n ,(,) = 0, (3.47) 

where D (z) is the straightforward extension of Dq(t), and 

C (z):=\^\ = r positive, 

[ Co{T) tor z = — t negative. 

Imposing the Frobenius expansion, eq. (3.33), as the starting point for the numerical evo- 
lution, which correponds to the boundary condition ^f n (z = +00) = 0, we were able to nu- 
merically solve the equations of motion (3.47), evolving them backwards to z — > —00. Note 
that the coefficient functions and ttt\ are continuous, while their first derivatives have 

D (z) D (z) ' 

a discontinuity at the origin 7 . In accordance with general Sturm-Liouville theory, we were 
able to find a discrete spectrum of (axial-) vector mesons by imposing ^ n (z — > —00) = at 
the opposite boundary. 

Numerical results. There will be two distinct cases according to the large z asymptotics, 
i.e., depending on whether the asymptotic wave functions approach zero from above or below: 

• Vector mesons (V): Two identical types of solutions for lim^i^oo ty 2 n-i(z) = + or CP. 

• Axial- vector mesons (A): Two identical types of solutions corresponding to 
lim z ^ ±00 # 2n (z) = ± or 0*. 

The corresponding eigenvalues A n are listed in table 1. Note that, in accordance with the 
discussion above, there are two types of solutions for small \z\, namely the "cosine-like" and 
the "sine-like" solutions and there is an approximate parity symmetry. Again, the deviations 
from parity invariance are due to Aj(t) 7^ A 7 (t). 

We found solutions *$? n (z) and the corresponding eigenvalues \ n for the first 20 eigenvalues. 

The results for a selection of wave functions are shown in figure 2. 

There is also a non-renormalizable zero mode subject to d z ^o(z) oc D(z)~ l and 

lim^ioo ^o(z) = ±1, associated with the pions in the theory (see discussion in [40], sect. 

4). This mode is shown in figure 3. 

It is possible to fix the mass scale M* by identifying the lightest vector meson with the 
p-meson p(770) of mass M p = 775MeV, yielding Ml = = (2141MeV) 2 . Similarly, the 



7 In physical applications, a discontinuity usually signals an abrupt change in the propagation medium of 
a wave and can be treated be demanding appropriate "junction conditions" across the discontinuity. 
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n 


1 


2 


3 


4 


5 


6 


A 2 n-1 - M 2 


0.1310 


0.4785 


1.1081 


2.0267 


3.2238 


4.6948 


A - M 2n 
^2n - -JJT 


0.2582 


0.7785 


1.5552 


2.6004 


3.9217 


5.5214 



Table 1: Some numerical values for the dimensionless (axial-) vector meson masses. 

lightest axial vector meson can be identified with ai(1260), the second lightest vector 
meson with p(1450), and so forth. Moreover, it would be possible to fix k so as to 
obtain a realistic value for the pion decay constant. In table 2, we compare the mass ratios 
of the first few mesons in our model (DKS) (for the antipodal embedding tq = 0) with the 
corresponding ratios found in the Kuperstein-Sonnenschein (KS) [40] and Sakai-Sugimoto 
(SS) models [35] and experiments [48]. 



k 


( Afe+i \ 

V Al /dks 


V Al Jks 


( Afe+i \ 

\ Al /ss 


\ / cxp. 


1 


1.97 


2.68 


2.34 


~ 2.51 


2 


3.65 


5.36 


4.29 


~ 3.56 


3 


5.94 


8.88 


6.79 


~ 4.54 


4 


8.46 


13.3 


9.85 


- 4.15 



Table 2: Ratios M% +1 /Mf = A^+i/Ai for k = 1, . . . , 5. The numerical values in the DKS case 
are obtained for the antipodal configuration r = 0. One should note that the DKS ratios 
are closer to the experimental values than other holographic models, while for k > 3 neither 
model matches the experiments well. 
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Figure 2: (Normalized) wave functions y/K^ n {z) for the first three vector mesons n = 1, 3, 5 
and the first three axial- vector mesons n = 2,4, 6. 



3.2. Scalar mesons 

In this section, we are interested in studying the scalar fluctuations around the flavor brane 
embeddings. As for the vector mesons, we will derive a five- dimensional effective action 
up to quadratic order for the 5d scalar fluctuation modes corresponding to (pseudo-) scalar 
mesons. The potentially relevant modes stem from the Kaluza-Klein reduction of the 8d 
scalars fluctuations 59\, 5<pi transverse to the /J7-brane embedding and scalar modes of the 
8d gauge field, i.e., SA 5 , SA 6 , 5A 7 , respectively. 

It is important to note 8 that, for the antipodal configuration tq = 0, there is an unbroken 
£7(1) that rotates the S 2 around the two antipodal points, but leaves the JJ7-branes invariant. 
In this extremal case there is also an emergent U(1)r symmetry that shifts the angle ijj by 
a constant, acting on the 7J7-branes but leaving the induced metric invariant. While the 
unbroken U(l) symmetry may rotate fluctuations of $i into fluctuations of 0i, we will see 

8 We would like to thank A. Dymarsky for clarifying this point to us. 
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Figure 3: Zero 


mode a/k\I/o(- 2; )- 





that from the expansion of the DBI-WZ action there are no mixing terms between the 89\ 
and S(pi fluctuations and, interestingly , the scalar fluctuation d~A 7 . This happens accidently 
only for the antipodal case, where <j)'\{r) = classically 

The relevant fluctuations for a QCD-like theory are singlet states 9 under the 577(2) x 
U(1)r isometry of the internal S 3 . The simplest way to preserve the SU(2) x U(1)r isom- 
etry discussed above would be to disregard S6i, 5<pi, 5A 5 and SA 6 fluctuations, which would 
automatically break the isometry symmetry. It is possible that there exists a nontrivial 
singlet state, i.e., some linear combination of all the fluctuations. This will be left to a 
future investigation. In this subsection, we calculate the effective actions for all the scalar 
fluctuations. We will find non-trivial mixing terms between the 5<pi {o~9\) and the (5Aq) 
fluctuations 10 . We will see that the 8A7 fluctuation does not mix with the other fluctuations 
suggesting its interpretation as a singlet and therefore a good candidate for a scalar meson 
in the dual theory. However, after performing a Kaluza-Klein expansion we do not find a 
consistent set of orthonormal eigenfunctions due to the discontinuity at the origin of one of 
the Sturm-Liouville coefficients arising from the DBI-WZ expansion. However, it is quite 
clear that this discontinuity is not a real physical effect (e.g., there is no charge located at 
the origin), but rather an artifact of our choice of coordinates. Therefore a smooth and 
stable spectrum of scalar mesons should exist. Unfortunately, we have not yet been able 
to numerically solve the corresponding Sturm-Liouville problem in a completely satisfactory 
way. 

3.2.1. Effective action for scalars fluctuations 
Consider the coordinate fluctuations 

<50i (z,t) = {2i:a')d 1 {x,T) , 

9 Non-singlet states with nontrivial charges under the SU(2) x U(1)r isometry have no counterpart in the 
dual QCD and will therefore be disregarded. 

10 Thcse terms arise from the metric components describing the non-trivial fibration of the S 3 over the S 2 
which were absent in the original treatment of Sakai and Sonncnschcin [39] . 
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>i(x,t) = (27ra')0i(x, r) . 



(3.49) 



As a consequence, the induced metric now contains linear and quadratic terms in the fluc- 
tuations 

£(5(8)) mi/ = Gg 1 g 1 d fl S9 1 d u 59 1 + GfafadftSfadvScf)! , 
fi(9(8))iMi = Ge^dfj.SOiSn + G^d^Scj)^ + SG^d^iSis , 

£(5(8) )ri = Go^drSOiSu + G^ l2 d T 6<f>iSi2 + £G0 l3 <9 T 5<Mi 3 , (3-50) 

where Gpg is the 10-d metric of the deformed conifold (see appendix B). On the other hand, 
the gauge field fluctuations 



8A 5 (x,r) 
SA 6 (x,t) 
SA 7 (x,t) 



(27ra')a5(x, r) 
(27ia')aQ(x, r) 
(2iTa')aj(x, r) 



(3.51) 



induce the field strength fluctuation 



£^23 
£^>i 



bd T 5(J) 



- + d T 8A 5 , 8F T 2 = d T 5A e , 5J> 3 = d T <L4 7 



coshr 
-<L4 5 , = SA 6 

d fl 5A 5 , ^2 = 6^5 A 6 , = ^5i 7 . 



5JT 12 = -5A 7 



(3.52) 



From the linear and quadratic terms in the DBI expansion we extract the actions for the 
scalar fluctuations 



S 



(6^,8 As) 
DBI 



f^7 

nil; 



d A xdrd 3 n^\E\\ h x l % iT (G^ - E^G^dnSfadJfa 



+ E L %5A 5 d u SA 5 - 2E L2 G (j)l2 d^ 1 dJA 5 

+ E TT [{G Ml - E 22 Gl i2 )(dMi) 2 + E u 

2E l2 G, l2 dMi (- b -^ + d T 5A 5 
\ cosh T 

bd T 8d>\ 



bd T c 



coshr 



- + d T SA, 



2E 



t3 



E 12 5A, I - 



coshr 



d T SA 5 + E 22 G^ 2 SA 5 d T S<f>i 



+ E 22 E 33 (5A B ) 2 } 



(3.53) 



,5' 



(<50i,<5A 6 ) 
DBI 



J d A xdTd 3 n y ^\^h 1 / 2 7 1 ^[(G eiei -E u G 2 ei ^5e 1 dje 1 



E 2 %5A 6 d v 5A 6 + 2E 12 G Bll d ll 50 l dJA 



E 1 



(G 



0101 



# Gfti)(0rtf0i) + E Z2 (d T 5A 6 ) 2 + 2E L2 G ei id T 5A 6 d T 56 



- 2E t3 [E 12 SA 6 d T SA G - E n Ge 1 i5A 6 d T 59 1 ] 
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£ n £ 33 (cL4 6 ) 2 }, (3.54) 



S^f = -0 J d A xdrd 3 n^/\E\^ - 2E r3 d T 5A 7 + 2E 12 5A 7 
+ E 33 h 1/2 ^ v d^A 7 d u 5A 7 + E TT E 33 (d T 8A 7 ) 2 

+ E U E 22 (5A 7 ) 2 -2E 12 E r3 (d T 5A 7 )8A 7 y (3.55) 

The tensor E MN is defined as the inverse of Emn = (g(g))MN + J~mn (see appendix B for 
details). The WZ action is given by 

Swz = -y J ^A^AC (4) 

= — —^-tappa i h~ 1 J r a/3 J r pa d 4 xdTdQ 1 d 3 Q (3.56) 
8ds J 

where a = (r,fl l ). This action contributes to the scalar fluctuations as 
s Va 5 ,sa 6 ,sav) = -tl J d^xdrd^nh-^ - J r rs5A 7 + 7 12 d T bA 7 

- 6A 5 (JA^l + 9t sa 5 ) - (d T SA e )SA 6 - (d T 6A 7 )6A 7 \ , (3.57) 
\ cosnr / J 

As expected, the linear term arising from the WZ action cancels with the linear term 
arising from the DBI action. The results for the quadratic terms show that there is a non- 
trivial mixing between the scalars <5</>i (S9i) and 5A 5 (5Aq). A linear combination of these 
fluctuations is not enough to diagonalize the action. A simple solution to this problem is to 
set 5<fri = 89i = 5A 5 = 5A e = but we don't discard the possibility of a non-trivial singlet 
combination. 

3.2.2. Scalar fluctuations of A 7 

The only remaining scalar mode of interest is the one associated with fluctuations of A 7 . 
Summing the DBI and WZ terms we can extract the action for the fluctuation a 7 

Sties. = -*J dx 4 dr r-G{r) rT <Vr(^, r) d u a 7 (x^ r) + M 2 H(r) d T a 7 (x», r) d T a 7 (x^ r) 

+M 2 I(t) (a 7 (x», r)) 2 + M 2 J(r) (d T a 7 (x», r)) a 7 {x^ r)) , (3.58; 
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TT _ 1 



G(r) 
ff(r) 

J(r) 
J(r) 



^ 2 Vol S' 3 e 4//3 , M 2 = 6 rn ^7^2 as before and we have again utilized the 
unctions introc 

2K(t) coshr = 2D(t) 



where « = f/i 7 # s ^7ra, . , - 

rescaled dimensionless functions introduced in (3.12) to define 



h(r)K(r) 2 \ Al ~ h 



\ 



-i 



V 



coshr 
4 



9K 6 (r) cosh' 



3 if (r) cosh 3 r 



C{t) 



( 



h(r)K(r) 



A 7 



V 



coshr 3 if (r) cosh 3 r 
4 



x- 1 



h(r)K(r) cosh 3 r 



3if 3 (r) coshr 



where the second term in the first line of J(r) is a contribution from the WZ term . We will 
again drop the hats and superscripts for notational convenience and continue to use barred 
quantities H(t), I(t) and J(r) to denote functions defined on the left (anti-D-brane) branch, 
i.e., with the appropriate choice of A 7 substituted into the functions defined above. 
Strategy. The five- dimensional fluctuations on the two branches can then be decomposed 
as 

oo 
k=0 

The wave functions n (r) have to satisfy the following equations, 

k J drG(r)0 n (r)0 m (r) 

(G(r)r 1 \-8 T {H{r)d T ^ n {r)) + /(r)0 n (r) + J(r)<9 r n (r)] 

implying 



-A„0„(r). 



(3.59) 
(3.60) 



k J dr (H(r)d T (f) n (T) <9 T m (r) + /(r)0 n (r)0 m (r) + J(r) (<9 r n (r)) <f> m {r)) = ^\J r 



This will lead to a four-dimensional effective Lagrangian 



- / d 4 xJ2 ha 7 n \x^a 7 n) (xn + \ n Ml(a 7 (x»)) 2 



(3.61) 
(3.62) 



Therefore, a 7 represents a scalar meson with mass given by |A n M 2 . 

Again, we will study the D7-brane branch first and then generalize to the D7-brane branch. 
The equation of motion can be rewritten as 



H(t) H(t) J ' \2 n G(r) J H(t) 
22 



The leading large r behavior of the various coefficient functions is given by 
m = t lim m = 0, lim M = t lim J « - 



r^oo if(r) 3' T-+00 H(t) ' r^oo H(t) 9' t^co H(t) 

Therefore, it is easy to see that the normalizable solutions behave as 

n (r) = ±e~ 4+ s 5r <^(r), n (r -> oo) ->■ const. 
The resulting differential equation for n (r) reads 

J£(r) + S(r)&(r) + (T(r) + ^A n C/(r))^(r) = 0, 



4 
3' 



where 



S(t) 
T(r) 
U(t) 



H'{t) J(r) _ 8 + 4\/5 
#(r) ~~ H(t) 3 ' 

(4 + 2y/5) 2 _ 4 + 2v/5 / ff'(-r) J(r) 
9 3 \H(t) ~ W{t) 



H(t) 



T 



2Z(n 



(3.64) 

(3.65) 

(3.66) 

(3.67) 
(3.68) 
(3.69) 



Large r regime. It should be noted that at leading order in large r, the contribution of 
the background gauge field vanishes, since 

3 



lim (A 7 (r)-b(r)f 



lim d T (A 7 (r) - 6(r)) = 0. 



(3.70) 



The expansion of the relevant coefficient functions reads We can expand the coefficient 
functions S(t),T(t) and U{r) for large r as 



S(r) = J2s k (r)e- 2k \ T(r) = J>(r)e- 2fcT , U(r) = e -§^ Mfc (r)e' 

k=0 k=0 k=0 

where the first few expansion coefficients are determined to be 
AV5 4 



2fcr 



(3.71) 



So(t) 
«i(t) 
to(r) 
ti(r) 
Uq(t) 



5 + 4r' 



16 + - i 



4 + 2^5 4 

3 5 + 4r' 
4 + 2^ 



16(1 



32 /i 



2 5/3 f 



15 3 



Ml (r)=2 5 / 3 



V 8 2 
9 9 
16 ~ 2 



250 



(-96 + 555r + 300r 2 ) 



2 8/3 ^_ 



15 



r (l-2r) 
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Here we included at least the leading order in r. Expanding 



, _2fcT 

Tie 3 , 



(3.72) 



fc=0 



and using (3.65), it is easy to find a recursion relation for f3 n ^ 
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m=0 



fc-i l 



+ 77 a w M m /? ; 



n,fc-3m-l 



0. 



m=0 



Setting /3 nj0 = 1, we find the first non-trivial coefficient to be 

= 27 2~ 13 / 3 A„ r (71 - 31\/5 + 2(-2 + >/5)r) 



(3.73) 



(3.74) 



Numerical results. In principle, one could proceed as before and try to solve the Sturm- 
Liouville problem via a shooting method. However, we did not obtain a stable numerical 
solution for the eigensystem of eq. (3.63). This is due to the discontinuity of the coefficient 
S(t) at the origin t = \z\ = (related to our choice of coordinates). We tried to reformu- 
late the problem slightly, using two different eigenvalues for the two branches and glueing 
the solutions smoothly at the origin. However, in that case the set of solutions do not sat- 
isfy orthonormality. Therefore we were not able to find a stable spectrum of scalar meson 
fluctuations for a-j on the full domain of the flavor D7-D7 branes. 



3.2.3. Scalar fluctuations of 61 (and (j>i) when 5A§ = SAq = 0. 

Finally, we would like to analyze the effective action for the scalar fluctuation 8\ in the 
special case when the gauge field fluctuations are turned off, i.e. 5A 5 = and 5A 6 = 0, 
which is consistent with the (A)SD eqs. (2.23) only for the case r = 0. For this case, the 
effective action associated with fluctuations around 8^ 



I takes the form 



where 



<7 

j 1 



5d,eS. 



d A xdr {A^ifd^duSOi + M 2 5(r)(d r 50i) 2 } 



A(t) 
B(t) 



(±7 

h(r) cosh 2 t 



Vo\(S 3 )e^(g s Ma') 2 
2 



3K 2 ( 



(A 7 - by - Ah(r) = C(t)D(t) - Ah(r) 



K 2 (t) cosh 2 r 



D 2 (r) 



Mr 



1 - Ah(r) h(r) cosh 2 r + 



D(r) 



3K 2 (t) 



(a 7 - by 



C(r) 



(3.75) 



(3.76) 
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on the D7-brane. On the D7-brane branch, one has to replace Af{r) by A-j{r). We have 
made use of the same dimensionless functions defined in (3. 12). The effective action for the 
5<f)i fluctuations yields a similar result. The five dimensional fluctuation 6'i(x M ,r) can be 
decomposed as 



(3.77) 



n=l 



so that we get a four dimensional effective action 



J 4d,eS. 



n=l 



4£ / d ' x rrd^W^ + M^Xniei 



(")\2 



(3.78) 



where the normalization conditions and equations of motion for the wave functions £ n (r) 
read 



k J dr A(t)U(t)Ut) 
iA^y'dAB^drUr)) 



AnCn(T) 



(3.79) 

(3.80) 



with X n := M%/M% denoting the eigenvalues of the scalar fluctuations. Again, we are solving 
a Sturm-Liouville type problem and in principle we can follow the same strategy as for the 
vector mesons and use a shooting method to solve the equations of motion, rewritten as 



Br) 



Br) 



(3.81) 



where 



B'(r) 
B{r) 

A(r) 



C(r)D'(r) - 2h(r) 

C(r)D(r)-4h(r) 
C(r) h(r) , 2 (A 7 (r) - b{r)f 



h'( 



D'( 



C'i 



h(-, 



D( 



C(r) 



B(r) D(j) K 2 (r) 3 K 4 (r)cosh 2 r 



(,3.82a) 
(3.82b) 



on the _D7-brane branch. Again, on the D7-brane branch, one has to replace A 7 {r) by A 7 {r). 

Thus we again find a discontinuity at the origin r = \z\ = , this time of the coeffi- 
cient B'(t)/B(t). As a consequence we did not obtain a stable numerical solution for the 
eigensystem of eq. (3.81). 



4. Conclusions and Outlook 

In this paper, we studied the vector and scalar mesons associated with the vector and scalar 
modes of the type IIB supergravity solution, presented in [33], corresponding to a stable 
non-supersymmetric D7/-D7-brane embedding geometrically realizing flavor chiral symmetry 
breaking in the Klebanov-Strassler background. In [39], it proved difficult to find a solution 
to the equations governing the background gauge field. Therefore the authors set A a = 
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0, a = 5, 6, 7, even though this is not allowed since the background gauge field couples to 
a non-trivial current. A resolution of this problem was presented in [33] (see discussion 
above), enabling us to study all relevant vector and scalar fluctuation modes (at least for 
the antipodal embedding tq = 0), including the pseudo scalar mode associated with the 
non-trivial background gauge field A?(t). 

We found a consistent set of eigenfunctions for the vector mesons, axial vector mesons 
and the pion. From the corresponding eigenvalues we estimated the masses and found a 
reasonable agreement with previous models and with experiment. For the scalar fluctuations 
we calculated the effective action from the DBI-WZ expansion but didn't find a singlet 
combination that leads to a consistent set of eigenfunctions. 

The meson wave functions that we obtain in this model are only approximately parity 
symmetric in the radial coordinate z = ±r. The parity symmetry in the ^-direction is 
broken, due to the presence of a non-trivial background gauge field that turns out to be 
asymmetric in the ^-coordinate. This issue deserves some further elaboration: 
The charge C and parity P symmetries of the five- dimensional fields are inherited from the 
7+1 dimensional DBI action for the D7-branes. After Kaluza-Klein reduction on S 3 , they 
become the five- dimensional symmetries P : (x M , z) — > (— x M , — z) and C : A^{x^,z) — > 
— A^(— x M , —z). The classical configuration that minimizes the 7+1 dimensional action (i.e., 
the .D7/-D7-embedding together with an ASD/SD gauge field on its world volume) breaks 
the z — > —z parity symmetry. However, it is important to observe that parity symmetry 
is restored asymptotically at large z. Therefore we can define the parity quantum numbers 
of the vector and axial-vector mesons according to their coupling to external gauge fields 
V = A L + A R (representing a photon) and A = A L — A R , resp., at the boundary \z\ — > oo. 
It seems that the asymmetry of the wave functions in the ^-direction could be related to a 
violation of charge conjugation in 3+1 dimensions, but we leave the precise physical inter- 
pretation to future research. 

It would be very interesting to study the meson spectra for general tq > in order to 
be able to study more realistic particle phenomenology. However it seems hard to solve the 
field equations for the background gauge fields even numerically in that case. We leave this 
as an open problem for future study. Some of the present authors are currently investigating 
baryons in this setup, which is also relevant phenomenologically, e.g., for studying nuclear 
physics in the context of this holographic model. The baryons in this context are usually 
modelled as instantons, where the instanton number is interpreted as the baryon number. 
It seems promising that, unlike in the Sakai-Sugimoto model, the baryon size does not go to 
zero as A -1 / 2 in the large t'Hooft coupling limit A — > oo, which is related to the fact that the 
five-dimensional effective k does not depend on A. It would also be of considerable interest 
to study in detail other phenomenological questions such as pion and vector meson form 
factors and deep inelastic scattering of pion and mesons in the context of this model. 
Another important aspect in the context of holographic QCD is the introduction of holo- 
graphic baryons to study, e.g., nucleons and nuclear forces in holographic models. For a 
more recent discussion of holographic baryons in the DKS model, cf. [38]. This line of 
work is based on original ideas by E. Witten [49] and has been pursued by many authors in 
the Sakai-Sugimoto D4/D8 model [50,51,52,53]. The conclusion of a recent article [54] was 
that there cannot be a 'realistic' combination of attractive and repulsive nuclear forces in 
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the Sakai-Sugimoto model. As explained in [54], in order to get an attractive nuclear force 
slightly greater than the repulsive nuclear force, 

V iep ir) oc +r ~ 1 e ~ rM ^-, \/ attr -(r) oc -r^V^-. (4.1) 

one needs the lightest isoscalar scalar meson to be lighter than the lightest vector meson. 
Unfortunately, we have not managed to find a stable scalar meson spectrum in our setup for 
To = 0. The situation may be different at tq ^ and especially in the limit To — > oo where 
we expect the relevant equations to simplify significantly. Therefore, we can not yet directly 
address the question of the mass of the a-meson compared to the cu-meson. 

In conclusion, it would be worthwhile to study in detail the non-antipodal embedding 
in the DKS model (for To > 0), to see if it can provide more realistic physics. While some 
phenomenological properties of the DKS model look very promising, more work must be 
done to see if the model can address phenomenological questions in QCD successfully. 
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A. The deformed conifold 

The deformed conifold is a non-singular deformation of the singular conifold: 

4 

z i = -ZdeUjW = e 2 , (A.l) 

i=l 

where 

W = ( Wu Wl2 \ = — ( Zs + iZ4 Zl ~ iZ2 ^ (A 2) 

In terms of the radius p of the S 3 , TrW^W = p 2 , one can define a new coordinate, 

p 2 = e 2 coshr. (A.3) 

This is useful, because r = at the tip of the deformed conifold, where the S 2 shrinks to 
zero. There are many different coordinate bases used to parametrize the deformed conifold 
in the literature. The reader is referred to the appendices of the review article [27] for a 
detailled discussion of the different coordinate systems, including some of their shortcomings. 
Klebanov and Strassler [18] use the diagonal basis originally introduced by [56,57]. The 
metric reads 

^ 4/3 d4) = (A.4) 

^ (s^W^ 2 + {g5)2) + c ° sh2 © (( " 3)2 + {gr) + sinh2 © {{gl)2 + (/)2) ) ' 
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where 



K(t) 



(sinh(2r) - 2r) 



1/3 



(A.5) 



2 1 /3 s i n hr 

and the various R-R and NS-NS forms can be found, e.g., in [18]. For large r, in terms of 



the radial coordinate r 
singular conifold 



P 



2/3 



the deformed conifold metric approaches the metric of the 



(A.6) 



It was argued in [27], however, that this choice of vielbeins does not allow for the construction 
of a closed holomorphic three-form in the deformed conifold case, implying problems with 
the Calabi-Yau condition. Here, we will make use of the one-forms e, on S 2 and Ej on S* 3 
that were introduced by [55,28] and used by Dymarsky, Kuperstein and Sonnenschein [33] 
to write the metric (2.6). We utilize this basis in the main part of the paper. This basis is 
non-diagonal and one can observe that the S 3 is fibered over the S 2 explicitly: 



ei 

e 2 

El 
£3 



d0i, 

sin^id^i, 
cos ipc\9 2 + 



cos 9 



xacpx 



cub (fju.u 2 t sin ip sin 6 , 2 d0 2 = £1, 
— sin ipd92 + cos ip sin 6 l 2 d0 2 7 
d^> + cos 6^01 + cos6 l 2d02 = £3 

Since the non-BPS D7-brane embedding described in [33] covers the S 3 completely, i.e., 
preserves the SU(2)r invariance of the deformed conifold, it is not necessary for our purposes 
to use a basis in which the coordinates on S 3 are independent of the coordinates on S 2 . 
Such a trivialization basis was found in [58] and written down explicitly in [32] (see also 
[59]). It is stated here for completeness and because it may be useful in other contexts. In 



this S 3 x S* 2 basis, the metric becomes 



e - 4/3 d4) 



QK 2 (t 



idr 2 + h 2 3 ) + 



K{t) 



cosh ( — ) x 



h\ + hj + 4 tanh 2 




sm 




• (A.7) 



The one-forms hi 



1, 2, 3 manifest the fibration of the S 3 over the S 2 
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sin 9 
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(A.8) 
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The function K(r) is identical to the one defined above and the SU (2) left-invariant Maurer- 
Cartan forms Wi, % = 1, 2, 3 read 



— = — cos (5da + sin a cos a sin /3d/3 — sin 2 a sin 2 /M7, (A. 9) 

— = — sin /3 cos 7da — (sin a cos a cos /3 cos 7 — sin 2 a sin 7) d/3 

+ (sin 2 a sin (3 cos /3 cos 7 + sin a cos a sin /3 sin 7)d7, (A. 10) 

= sin (3 sin 7da + (sin a cos a cos (5 sin 7 + sin 2 a cos 7) d/3 

+ (sin a cos a sin (3 cos 7 — sin 2 a sin f3 cos /3 sin7)d7, (A. 11) 

where a,/?, 7 are coordinates on S* 3 . 

B. Expansion of the DBI action 

Consider the pullback of the 10-d coordinates X = (x^, r, 9±, $1, fij) on the D7 brane 
worldvolume with coordinates X a = (x^,r, Qi) where 8\ = #i(r, x), 0i = 0i(r, x) and 
dQi = €i (i = 1,2,3) 

The pullback is defined by 

(<?( 8 )W = G PQ d M X p d N X Q (B.l) 
where Gpq is the 10-d metric defined by 

A/3 



cfe 2 = h~ l,2 r]^dx^dx v + /i 1/2 — A cosh r|d# 2 + sin 2 M</>? + dtt\ + rffi 2 

+ — ^— (dMOa + Bin0i#idn 2 ) + or „ 2 : ^r 2 + (rffi 3 + cosM0i) 2 ]j (B.2) 
cosh r 3A 13 cosh r J 

For the antipodal case the pullback takes the form 

(g(8))MN = (h- 1 /\ v ,h 1 / 2 g a p), (B.3) 

where a = (r, 1, 2, 3) and 

e 4/3 / 2 2 \ 

9af> = — Kc0shT diag 1 3^3coshr ' *' 3^ cosh r J " (R4) 



On the other hand, the gauge-invariant field strength is 

d T A 5 ) dr A ei + d T A 6 dr A e 2 + <9 r A 7 dr A e 3 



.F = P[B {2) ] + F 
bd T 4>i 



cosh r 

- A 5 e 2 Ae 3 + A 6 e 1 Ae 3 + (b- A 7 )ei A e 2 + 

+ <9 M A efa^ A ei + O^Aq dx^ A e 2 + <9 M A 7 dx^ A e 3 + -F^dx 11 A rfx 1 ' 

=: ^J^mn dX M A dX N (B.5) 
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where 2-ira' has been absorbed in the gauge fields. In the antipodal case it takes the form 

Fmn = (0,-Fa/s) , (B.6) 
with ^ 

J^12 — b — A 7 , T T \ — J>2 = T\z = -^13 — , J>3 = J 7 12 , (B-7) 

The DBI action can be written as 

Sdbi = -ih J d 4 xdrd 3 tle-*y/\ det(E MN + 8E MN ) \ . (B.8) 

where 

Emn = {g(8))MN + Fmn , (B-9) 

is the background matrix and 

SE M N = <%(8))a/AT + ^AfJV (B.10) 

the fluctuation matrix. The background matrix has an inverse given by 

E M N = (h 1 / V, E a ?) (B.ll) 

where 

E «p = h ~ 1/2 V^~\ g -i _ h -i/2 g -ij: g -i}c$ (B . 12 ) 
\E\ 



and the square root of its determinant can be written as 



E\ = VW\ + h- 1 = l^^h- 1 [7\ 2 + g ll922 h) (B.13) 

V 911922 

We omit for a moment the tensor notation. We want to expand the determinant 

det(E + 8E) = det(E) det(l + E~ X 5E) , (B.14) 

Using the identity 

det(l + X) = exp{tr[ln(l + X)]} (B.15) 

and the logarithm expansion 

ln(l + X) = X-^X 2 + ..., (B.16) 

we obtain 

det{E + 5E) = det(£){l + trfE" 1 ^] + ^(trfE" 1 ^]) 2 

- ^trKE^Efj + OdSE) 3 )}, (B.17) 

so that 

y/-det(E + 6E) = ^/^det(E)\l + hT[E- 1 5E] + l(tT[E~ 1 5E}) 2 

L 2 8 

- ^[{E^SEf] + 0((5E) 3 )j . (B.18) 
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